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ABSTRACT 



Context. Understanding the relation between the star formation rate (SFR) and stellar mass (M+ ) of galaxies, and the evolution of the 
specific star formation rate (sSFR=SFR/M lr ) . 

Aims. We examine the dependence of derived physical parameters of distant Lyman break galaxies (LBGs) on the assumed star for- 
mation histories (SFHs), their implications on the SFR-mass relation, and we propose observational tests to better constrain these 
quantities. 

Methods. We use our SED-fitting tool including nebular emission to analyze a large sample of LBGs from redshift z ~ 3 to 6, assum- 
ing five different star formation histories, extending our first analysis of this sample (de Barros et al. 2012, paper I). In addition we 
predict the IR luminosities consistently with the SED fits. 

Results. Compared to "standard" SED fits assuming constant SFR and neglecting nebular lines, assuming variable SFHs yield sys- 
tematically lower stellar masses, higher extinction, higher SFR, higher IR luminosities, and a wider range of equivalent widths for 
optical emission lines. Exponentially declining and delayed SFHs yield basically identical results and generally fit best. Exponentially 
rising SFHs yield similar masses, but somewhat higher extinction than exponentially declining ones. We find significant deviations 
between the derived SFR and IR luminosity from the commonly used SFR(IR) or SFR(IR+UV) calibration, due to differences in the 
SFHs and ages. Models with variable SFHs, favored statistically, yield generally a large scatter in the SFR-mass relation. We show 
the dependence of this scatter on assumptions of the SFH, the introduction of an age prior, and on the extinction law. We show that 
the true scatter in the SFR-mass relation can be significantly larger than inferred using SFR(UV) and/or SFR(IR), if the true star 
formation histories are variable and relatively young populations are present. Different SFHs, and hence differences in the derived 
SFR-mass relation and in the specific star formation rates, can be tested/constrained observationally with future IR observations with 
ALMA. Measurement of emission lines, such as Ho- and [O n] /13727, can also provide useful constraints on the SED models, and 
hence test the predicted physical parameters. 

Conclusions. Our findings of a large scatter in the SFR-mass relation at high-z and an increase of the specific star formation rate 
above z ~ 3 (paper I) can be tested observationally. Consistent analysis of all observables including the rest-frame UV to IR and 
emission lines are required to establish more precisely true SFR values and the scatter in the SFR-mass relation. 

Key words. Galaxies: starburst - Galaxies: ISM - Galaxies: high-redshift - Ultraviolet: galaxies - Radiative transfer 



1. Introduction 

One of the important, recent results of multi-wavelength galaxy 
surveys is the finding of a well-defined relation between the star 
formation rate (SFR) and stellar mass of galaxies - now often 
called the main sequence of star forming galaxies - and the 
evolution of this relation with r edshift, at least from the nearby 
Universe up to redshift z ~ 2 (jNoeske et a"fll2007l : iDaddi et"al] 
I2007t lElbaz et alj|2007t iRodighiero et alj|201 lb . The observed 
redshift evolution corresponds to an increase of the typical spe- 
cific star formation rate (sSFR) by more than an order of magni- 
tude from z ~ to 2, indicating a higher star formation activity 
of galaxies in the past. 

Recent analysis of Lyman break galaxies (LBGs), combin- 
ing often HST and Spitzer observations, have tried to examine 
whether a SFR-mass relation is already in place at higher red- 
shift, and how th e specific star formation rate of galaxies behaves 
at z ~ 4-7 (e.g. IStark et al.l l2Q09b ISchaerer & de Barros|l2010t 
lLabbe et alj l2010b iMcLure et alj |201 U IGonzalez et alj |201 lb 
The first studies have found a relation similar to the z ~ 2 



SFR-mass relation, which would indicate no evolut ion of the 
sSFR beyond redshift a plateau in sSFR - (e.g. lStark et al.l 
2009: IGonzalez et al]|201 lb . This result appears difficult to rec- 
oncile with most theoretical models, which successfully ex- 
plain the behavior of star formation properties from z ~ to 
2, but predict a continuing r ise of the specific s t ar formation 
rate towards higher redshift dBouche et alj l2010t [Dutton et al. 
2010: IWeinmann et alJl201 lblDave et alJl201 lb. although others 
propose alternatives (Krumholz & Deke]||2012l e.g.). More re- 
cent work, including different spectral energy distribution (SED) 
models, or allowing for a revision of dust attenuation have 
shown that the obse rved sSFR at z ~ 5-7 could be highe r than 
previously thought (Schaerer & de Barros 2010; de Barros et al 



2011 



2012). 



iBouwens et alJ 1201 lb lYabe et alJ 120091 Ide Barros et al 



If all galaxies at z ^ 2 obey a mass-SFR relation with star 
formation rate increasing with galaxy mass, it is evident that 
their star formation rate must increases with time. The appar- 
ent, small scatter in the SFR-mass relation, and the value of 
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the exponent a close to u nity in this relation SFR oc M? found 
by various studies ( e.g. iDaddi et alJ l2007t lElbaz et al] l2007t 
Gon zalez et al.l l201lb has led several authors to conclude that 
high redshift galaxies mu st have gone through a phase of quasi - 
exponential growth (e.g. iRenzinil 2009; Mai-astojietalj|2010). 
Earlier studies have independently advocated rising star for- 
mation histories (SFHs) for high redshift galaxie s from hydro- 
dynamic simulations and semi-analytical models dFinlator et al.l 
120071 120 10t iLee et alj|2010t iFinlator et al1l201 lh . Arguments in 
fa vor of rising SFHs, at least on average, haye bee n put forward 
bv lFinkelstein et al.l(l2010h : lPapovich et all (1201 lh to explain the 
evolution of the UV luminosity function. If representative of the 
star formation history of individual galaxies, such star formation 
histories would indeed by quite different from exponentially de- 
clining or constant SFHs, most frequently assumed in the litera- 
ture to model/analyze the observed SED of distant galaxies. 

Assumptions on the star formation history affect the phys- 
ical parameters of galaxies derived from SED fits. This 
is therefore gener a lly treated as a free parameter. See e.g . 
Red dv et al.l d2010l) : lMaraston et al.ld20Tol) ; lWuvts etaTJ (1201 lh : 
Redd v et all (1201 21 for examples on z ^ 2 star forming galax- 
ies. Determinining star formation histories and their associated 
timescales is also important, as this may provide constraints on 
different modes of star formation and on feedback processes at 
high redshift, recognized as key features for our understanding 
of galaxy evolution. For example, star formation may proceed on 
different timescales if regulated by cold-accretio n or by star for- 
mation feedback or triggere d by mergers (cf. iKhochfar & Silk! 
1201 U IWvifhe & LoeblbOl lh . These points already illustrate the 
interest and importance of determining the star formation histo- 
ries and typical timescales of star formation at high z. We here 
wish to discuss and reexamine these issues on the basis of an 
analysis of a large sample of LBGs covering redshifts from ~ 3 
to 6, and to present observational tests which should be able to 
distinguish different star formation histories. 

One of main arguments often invoked to argue for rising 
star formation histories is the small scatter of the SFR-mass 
relation. However, it should be recognized that generally the 
SFR is derived from observables - the UV or IR luminosity 
- which depend on relatively long 100 Myr) timescales, 
and which are assumed to be at an equilibrium value, which 
is only reached after this timescale and for constant SFR. With 
these assumptions entering e.g. the common ly used SFR(UV) 
or SFR(IR) calibrations or iKennicutl dl998h . it is natural that 
the "observational" scatter is smaller than the true scatter in the 
current SFR, if typical ages are less than 100 Myr and/or the 
timescale shorter than this. Especially at high redshift, where 
timescales are shor ter (e.g. the dynamica l timescale decreases 
with (1 + zT 3/2 , cf. Wyithe & Loebll20Tll) one should therefore 
carefully (re)examine the SFR-mass relation and its "tightness" 
using consistent diagnostics, and with all the available obser- 
vational constraints including on age, star formation timescales 
and histories. In any case, the conclusion that rising SFHs are 
favored is a priori inconsistent with the assumptions leading 
to it. Although this has been re-examined for z ~ 2 sam ples 
dMarastonetalJ2010t IWuvts et alJioTTt iReddv et al.ll2012b this 
obviously calls for a revision at higher redshifts. 

Most SED studies of LBGs at z ^ 3 have assumed con- 
stant star formation rates, or exponentially declining SFHs 
to determine the p hysical parameters of these galaxie s (e.g . 
Egamiet al. l [20051: ISchaerer & Pellol 120051: lEyles et al] [2 005 ; 
Verma et al.ll2007l: lYabe et al.112009): IStark et al]l2009t ILee et al j 
20 10t ISchaerer & de Barrosl 120101: bonzalez et all 1201 lb . some 
of them imposing no dust attenuation, motivated by the blue UV 



sl opes observed at high redshift. Notable examples are the work 
of Finlato r~et al.l (l2007h who analyzed 6 galaxies at z > 5.5 with 
different SFHs, including rising ones taken from they hydro- 
dynamic simulations. The physical parameters they derive are 
consistent with those using simple parametrised histories, albeit 
with reduced uncertainties. Their study also shows that SED fits 
with rising star formation histories yield a higher attenuation 
than inferred assuming constant S FR, a result no t yet a ppreci- 
ated enough, which we also find in lde Barros et al. (2012) and in 
this paper. 

Another drawback of most earlier studies is that the contri- 
bution of nebular emission (most emission lines) is not taken 
into account, an effect which can significantly alter the ages, 
masses, and other physical parameters derived from SED fits, 
as demonstrated by Scha erer & de Barrosl (120091 l20l6T) . Indeed, 
there is now clear evidence for the presence of nebular lines af- 
fecting the broad-band photometry of high-z star-forming galax- 
ies (Ly man alpha emitters and LBGs in particular), as discussed 
e.g. in ISchaerer & de Barrosl (1201 lh . The best demonstration 
comes from LBGs with spectroscopic redshifts between 3.8 and 
5, among which a large fraction shows a clear excess in the 3.6 
/mi filter with respec t to neig hboring filters (K and 4.5 fim), as 
shown by Shim et al. I J201 lab . At these redshifts Ha is located in 
the 3.6 /mi filter, whereas very few lines are expected at 4.5 /im, 
showing that the observed 3.6 /im excess is naturally explained 
by strong Ha emission. 

Given these limitations of published SED studies and the 
interesting results o btained from our work on a small sample 
of z ~ 6-8 LBGs (ISchaerer & de Barrosl 1201 oh . we have re- 
cently undertaken an extensive study of the physical param- 
eters of a large sample of ~ 1400 LBGs at z ~ 3-6, us- 
ing our state-of-the art photometric redshift and SED fitting 
model including nebular emission, and considering a range of 
different star formation histories. Among the numerou s de- 
tailed results obtained in this work dde Barros et al.ll2012l here- 
after dBSS12) we find in particular from the preferred mod- 
els: 1) a large scatter in the SFR-mass relation for the pre- 
ferred, variable star formation histories, 2) higher dust atten- 
uation than obtained from models assuming SFR=const and 
from standard methods using the UV slope, and 3) a higher 
sSFR than commonly obtained, and a rising sSFR with redshift. 
Our models therefore reconcile the observationally-inferred spe- 
cific star formation rate with predictions from cosmological 
models predicting a continuous rise of the sSFR with redshift 
dBouche et al.ll2010t iDutton et al.ll2010t IWeinmann et al.ll201 it 
iDave et al.l 1201 it iKrumholz & Dekell 12012b . The results from 
dBSS12 have also other important implications, e.g. on the typi- 
cal star formation timescales at high redshift, and on the cosmic 
star formation rate density. 

Given the importance of these findings, it is of interest to 
carry out additional independent tests of the models and to pro- 
vide further constrains on the ages, extinction, and star formation 
histories of high redshift galaxies. In this paper we present pre- 
dictions allowing such tests. Furthermore, we extend the study 
of dBSS12 by exploring other star formation histories, not con- 
sidered in our previous paper. I ndeed, while d BSS 1 2 adopted 
the fixed, rising star history of Finla tor et al.l {2011) obtained 
from the hydrodynamic simulations, we here explore exponen- 
tially rising SFHs and so-called "delayed" histories with variable 
timescales. These histories were previously app lied to the SEP 
fits of galaxies at lower redshift, e.g. at z ~ 2 by Maraston et al.l 
(l2010b : IWuvts et al.l d2011b : lReddv et all d2012b and found to be 
preferred over other simple star formation histories. 
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First, we examine the effect of the different SFHs on the de- 
rived physical parameters of LBGs. We then present the implica- 
tions these different model assumptions have on the SFR-mass 
relation. After that we present consistent predictions for the in- 
frared luminosity of these galaxies, based on their SED fits and 
the assumption of energy conservation, abandoning the standard 
SFR-Lir conversions. We show in particular that such a consis- 
tent prediction yields a smaller scatter in the observables (Lir) 
than in the current SFR. The true star formation rate of LBGs at 
Z ^ 3 could thus very well show a large scatter around a "main 
sequence" even if the UV and/or IR luminosities show a small 
scatter. Our models also show that star formation histories can, 
to some extend, be distinguished from measurement of IR lumi- 
nosities, which will become possible with ALMA, since differ- 
ent amounts of UV attenuation are expected for different SFHs. 
We also present the predicted strength of some selected emission 
lines, which can be used to test our models and the different star 
formation histories. 

Our paper is structured as follows. The observational data 
and the method used for SED modelling are described in Sect. 
[2] The dependence of the physical parameters on the assumed 
star formation histories is discussed in Sect. [3] Our general pre- 
dictions for the IR emission are presented in Sect. [4] Specific 
predictions for z ~ 4 - 6 LBGs and ALMA are given in Sect. 
[5] The predicted strengths of optical emission lines are shown in 
Sect. [6] We discuss our results in Sect. [7] and Sect.[8]summarises 
our main conclusions. We adopt a A-CDM cosmological model 
with Z/ =70 km s _1 Mpc~', Q„,=0.3 and Q A =0.7. 



2. Observational data and SED modelling 



Ide Barros et alJ d2012l hereafter dBSS12) have analysed a large 
sample of z ~ 3-6 dropout-selected galaxies in depth using an 
up-to-date photometric redshift and SED-fitting tool, that treats 
the effects of nebular emission on the SEDs of galaxies. In their 
homogeneous analysis they determine the main physical proper- 
ties, such as the star formation rate (SFR), stellar mass, age, and 
reddening. They assess carefully their uncertainties, and discuss 
the evolution of these properties with redshift. We here extend 
these simulations to include other star formation histories and 
we examine the predicted IR luminosities and equivalent widths 
of selected emission lines from these galaxies. 

2.1. Photometric data and sample selection 



We ha ve used the GOODS-MUSIC catalogue of ISantini et ail 
(2009), which provide photometry in the U, B435, V606, 1776, 
zssolp, J, H, and K bands mostly from the VLT and HST, and the 
3.6, 4.5, 5.8, and 8.0 fim bands from th e IRAC cam e ra onb oard 
Spitzer. Us i ng s tandard criteria as in iNonino et al. (2009) and 
Stark et all (120091) we then selected U, B, V, and i-drop galaxies. 
To reduce the contamination rate (typically ~ 10-20 %) we only 
retained the objects whose median photometric redshifts agree 
with the targetted redshift range. This leaves us with a sample 
of 389, 705, 199, and 60 galaxies with median photometric red- 
shifts of z phot = 3.3, 3.9, 4.9, and 6.0. See dBSS12 for more 
details. 



modified to take nebular emission into account. In dBSS12 
we considered a large set of s pectral templates based o n the 
GALAXEV synthesis models of iBruzual & Charlotl d2003l) and 
covering different metallicities and a wide range of star forma- 
tion histories. A Salpeter IMF from 0.1 to 100 M was adopted. 
Our results can easily be rescaled to other IMFs with the same 
power law slope at high masses. Nebular emission from contin- 
uum processes and numerous emission lines were added to the 
spectra predicted from the G ALAXEV models as described in 
ISchaerer & de Barrosl ((2009), proportionally to the Lyman con- 
tinuum photon production. Following the results of dBSS12 we 
set the Lya line to zero for z ~ 3-5, and we assume the case 
B value for the i-drop sample. The intergalactic m edium (IGM) 
was treated with the prescription of M adaul ( 1 19951) . 

The free parameters of our SED fits are: redshift z, metal- 
licity Z (of stars and gas), the age f+ defined since the onset of 
star-formation, and a ttenuation Ay described by the Calzetti law 
(ICalzetti et aLl EzOOO). For comparison, a subsample of the data 
was also modeled using the SMC extinction law (cf. Sect. |7]i. 

Since one of the main objectives of this paper is to explore a 
variety of star formation histories (SFHs), we have adopted five 
different histories, which are summarized in Table [1] and illus- 
trated in Fig.Q] Note that for exponentially rising SFHs the SFR 
is set to zero, after a growth by more than 20 decades. The bolo- 
metric luminosity output per unit SFR, respectively its inverse 
SFR/Lboi, and the mass-to-light ratio in the V-band correspond- 
ing to these simple star formation histories are shown in Figs. [2] 
and [3] The behavior of these quantities are fundamental to ex- 
plain the main differences in star formation rate and stellar mass 
obtained from different models, discussed in Sect. [3] 

The first three cases, constant star formation, exponentially 
declining SF, and the a verage rising SFH p redicted by the hy- 
drodyn amical models of iFmlator et aTl (1201 ll) have already been 
used by Ide Barros et alJ (120121) . Two additional SFHs, exponen- 
tially rising, and so-called "delayed" SF showing a initial in- 
crease and subsequent decrease of the star formation, are also 
considered here. The way these assumed SFHs modify the de- 
rived physical parameters of LBGs and how they could be dis- 
tinguished is discussed below. 



Table 1. Star formation histories modeled in this work. All de- 
scribe the evolution of the SFR with time (age f). Three of 
them have one additional free parameter, the SF timescale t. For 
model B a minimum age of 50 Myr is assumed, for compari- 
son with other studies. Models A, B, and C are the same as in 
dBSS12. The color code used in the figures is indicated in the 
last column. 



ID 



SF history 



functional form parameters color 



A=DEC exp. declining 

B=REF constant 

C=RIS rising 

D exp. rising 

E delayed 



SFRoc exp^ 777 
SFR=const 
Finlator+2011 
SFRoc exp +,/r 
SFRcc fexp ' /T 



t, T 

t > 50 Myr 
t 

T, T 
t, T 



green 
black 
blue 
red 

yellow 



2.2. SED models 

Our SED-fitting tool, desc r ibed i n ISchaerer & de Barrosl {2009) 
and ISchaerer & de Barr os (2010), is based on a version of the 
Hyperz photometric redshift code of iBolzonella et ail ([2000), 



In practice we compute for each model set (i.e. star forma- 
tion history) SED fits for all combinations of z e [0,10] in steps 
of 0.1 Z = (0.02=Z o , 0.004, 0.001), t = (0.01, 0.03, 0.05, 0.07, 
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Fig. 1. Illustration of the star formation histories adopted in this 
paper. Shown is the SFR as a function of age. Within each model 
set A-E the curves are normalized to an arbitrary value. The 
colors used for the different SFHs are listed in Table Q] 

0.1, 0.3, 0.5, 0.7, 1., 3., co) Gyr for model AEl 51 age step s 
from to the age of the Universe (see iBolzonella et all uOOO). 
Ay e [0,2] mag in steps of 0.05. Minimisation over the entire 
parameter space yields the best-fit parameters and SED, along 
with other properties such as the stellar mass and star forma- 
tion rate (SFR) and UV magnitude. Other physical parameters, 
briefly described below, are determined from the best-fit SED. 
To determine the uncertainties of the physical parameters, we 
used Monte-Carlo simulations by running 1000 realisations of 
each object. From this we derived the probability distribution 
function for each parameter/quantity, either for each individual 
object or for (sub)samples. 

2.3. Predicted physical parameters and observables 

In addition to the fit parameters of our models, redshift z, metal- 
licity Z, age f, SF timescale r, stellar mass M*, and current star 
formation rate SFR, our models also allow us to determine other 
parameters and observables of interest. We briefly describe them 
here, for sake of clarity. Note that for all the physical parameters 
and observables we determine the detailed probability distribu- 
tion function (pdf) from the full MC simulations. From the ID 
pdf we derive in particular the median value and the 68% con- 
fidence interval, which we use to illustrate the uncertainties on 
these parameters. For most of the paper we use for simplicity the 
median value of the physical parameter of interest. 

2.3.1 . Emission lines 

From the numerous emission lines included in our spectral tem- 
plates we select some lines, primarily Lya, Ha, H/?, [O n] ,13727, 

1 For exponentially rising SFHs we explore r = (0.01, 0.03, 0.05, 
0.07, 0.1, 0.3, 0.5, 0.7, 1., 2., 3.) Gyr, for delayed histories r = (0.01, 
0.03, 0.05, 0.07, 0.1, 0.3, 0.5, 0.7, 1., 2., 3., 4., 5.) Gyr. 



Fig. 2. Predicted ratio SFR/Lboi as function of time for the some 
of the models with exponentially decreasing, delayed, and ex- 
ponentially rising SFHs (with timescales r = 10, 30, 300, 3000 
Myr indicated), compared to the case of constant SFR. The plots 
are shown here for BC2003 models with solar metallicity. Same 
color codes as in Fig.["D The standard SFR(IR) conversion factor 
from Kennicuttl (1 19981) is shown by the black dotted line. 



[O m] /L14959,5007, for which we save the emission line fluxes 
and equivalent widths predicted by the photometric SED fits. 
These lines can subsequently be used for comparison with spec- 
troscopic observations (existing or future). 

2.3.2. UV and IR emission 

For straightforward comparison with observations we use the ab- 
solute UV magnitude M1500 defined at 1500 A, and derived from 
the SED using a square filter of 300 A width centered on this 
wavelength. 

To define an emergent (observed) UV luminosity Luv we 
follow common practice by computing Luv = AFx> where 
A = 1800 A, and is the average flux between 1400 and 
2200 A. Assuming energy conservation, i.e. that the attenu- 
ated/extinct UV light is reemitted by dust in the IR, we can 
predict the amount of IR radiation expected from each galaxy. 
The corresponding IR luminosity Lir is simply computed from 
the difference between the intrinsic, unattenuated SED, and the 
best-fit SED with the parameters. In practice Lir is calculated 
from the attenuated energy integrated f rom 913 A to 3 fim (cf. 
Ida Cunha et all 120081; iNoll et all [2009). This quantity also in- 
cludes the nebular emission (lines and continuum) if present in 
the SED, consistently with our assumption of identical attenu- 
ation for both stars and gas. We have checked that the integral 
over the Balmer continuum and up to 3 //m yield basically iden- 
tical results. 

From the UV and IR luminosity just defined, we also com- 
pute Lir + uv, the emergent luminosity in the UV + IR domain, 
which is a good proxy of the total bolometric luminosity and to 
some extent also of the SFR (but cf. below). Similarly, we will 
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log(age) [yr] 



Fig. 3. Predicted mass-to-light ratio in the V-band, M+/Ly, for 
the same models as shown in Fig. [2] Same color codes as in Fig. 

m 

use the ratio of the IR to UV luminosity, Lm/Luv, a quantity 
kn own to correlate stro ngly with the amount of dust attenuation 
(cf. lMeureretai]|1999l) . 

2.3.3. Predicted IR, sub-mm, and mm fluxes 

Given the IR luminosity Lir predicted consistently from the 
amount of attenuation by dust, it is straightforward to predict 
the flux at various IR wavelengths, provided a corresponding 
spectral template for dust emission is adopted. We simply as- 
sumed standard modified black-body spectra described by F v oc 
v^ByiTd), with a fixed value of ft = 2 and three different dust 
temperatures Tj = (25,35,45) K, which span the range of ob- 
served dust temperatures in nearby and distant galaxies. The 
modified black-body is scaled to Lir computed from 8 to 1000 
fim (rest frame). From this, we determine the predicted IR fluxes 
in bands corresponding to observations with Herschel (PACS 
and SPIRE bands), several sub-mm bands (350, 450, 850, 870 
/mi), 1.2 and 2 mm bands, and ALMA bands 1-10 (correspond- 
ing to ~ 8.1 to 0.35 mm). 

3. Dependence of the physical parameters on the 
assumed star formation histories 

To illustrate and discuss the dependence of the physical param- 
eters on the assumed star formation histories we here choose the 
largest sample comprising 705 B-drop galaxies. The general be- 
havior discussed here is also found for most of the other samples 
at z ~ 3, 5, and 6. Some results for other redshifts are discussed 
in Sect. [5] 

3.1. Comparison of physical parameters 

It is well known that the physical parameters determined from 
broad-band SED fits depend on assumptions such as the star for- 



mation history. In dBSS12 we have shown how quantities such 
as the stellar mass, age, dust reddening, and SFR depend on this 
assumption and how they are modified with the inclusion of neb- 
ular emission. Here we briefly illustrate the dependence of the 
main physical parameters on SFH, including in particular the 
new histories explored in this paper (exponentially rising and 
delayed SFH; cf. TableQ]). 

Figures |4] to Q show how the M*, age, Ay, and SFR derived 
for different assumptions on SFH compare to models assuming 
exponentially declining star formation histories (our model A) 
for our largest sample, the 705 B-drop galaxies. This model is 
taken as a reference, since it generally provides the best fits for 
the vast majority of galaxies analyzed here (cf. below). For all 
models except for SFR=const (model B), we include the effects 
of nebular emission. Model B, called REF in dBSS12, serves as 
a "reference" model to ease comparisons with other studies in 
the literature assuming constant star formation, and neglecting 
nebular emission. When nebular emission is added to model B, 
this does not lead to very different parameters (cf. dBSS12). 



B-drops 




7 8 9 10 11 12 

log(M.) - decl [M ] 



Fig. 4. Comparison of the median stellar masses of our z ~ 4 
sample with 705 galaxies derived from SED fits assuming the 
star formation histories listed in Table Q] The masses are com- 
pared to those derived from model A, plotted on the x-axis. 
Black triangles show model B with constant SFR, blue and red 
symbols the rising histories (model C and D respectively), yel- 
low circles show the delayed SFHs (model E). For each model 
the same number of galaxies are plotted. Invisible data points are 
hidden close to the one-to-one relation. 

Figure|4]shows that most models yiel d similar stellar m asses, 
as already known from earlier studies (lYabe et al.ll2009i) . The 
largest differences are found with respect to models assuming 
constant star formation and a mimum age of 50 Myr (model 
B), which yield generally higher stellar masses. On average we 
find masses lower by a factor ~ 3 at z ~ 4 for declining SFHs 
(model A), where half of this effec t is due to inclusion of nebular 
emission (cf. lde Barros et all20 12). The difference is mostly due 
to the average shift towards younger ages, which have a lower 
mass-to-light ratio, as shown in Fig. [3] It is immediately appar- 
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Fig. 5. Same as Fig.|4]for the stellar age, defined since the onset 
of star formation. For each model the same number of galaxies 
are plotted. Invisible data points may be hidden, especially under 
the yellow symbols of model E plotted last, due to discrete age 
values. 
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Fig. 6. Same as Fig. [4] for the attenuation Ay. For each model 
the same number of galaxies are plotted. Invisible data points 
may be hidden, especially under the yellow symbols of model E 
plotted last, due to discrete age values. 



ent that the models with delayed star formation histories (model 
E) yields very similar results to model A. Exponentially rising 
SFHs with v ariable timescales (m odel D) and the rapidly rising 
history from lFinlator et all d201 lh (model C) give also compara- 
ble stellar masses, although with a larger scatter, when compared 
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Fig. 7. Same as Fig. |4]for the current star formation rate SFR. 



to model A. In short, all the models with variable SFHs and neb- 
ular emission considered here yield quite similar stellar masses, 
which are on average lower than those obtained from standard 
models without nebular emission. 

Age is notoriously known to be sensitive to m odels assump- 
tions (e.g. lPapovich et al.ll200lHYabe et al.l l2009). Indeed, abso- 
lute ages predicted by different models can differ by orders of 
magnitude, as shown in Fig. [5] Overall models with delayed SF 
(E) show again quite similar results to those with exponentially 
decreasing histories (A). The behavior of the rising SFHs (C and 
D) is quite comparable. By construction, model B yields only 
stellar ages > 50 Myr. Whereas for exponentially declining and 
delayed SFHs most objects are best fit with relatively short and 
comparable timescales t, the opposite is true for fits with expo- 
nentially rising histories. In this case the probability distribution 
function for t clearly peaks at the maximum value allowed here 
(r = 3 Gyr), which means that a relatively slow growth is pre- 
ferred for most galaxies. This seems quite natural given the find- 
ing that the rising SFH from Finlator et al ] (120111) fits on average 
less well than declining models. 

The attenuation predicted by the various models is shown 
in Fig. [6] As already discussed by dBSS12, constant SF models 
show generally lower dust attenuation than t he exponentially de- 
clining models, whereas the rising SFH of lFinlator et al.ld2011h 
predicts a higher attenuation on average. Rising star formation 
histories with variable timescales globally scatter around those 
assuming declining SF. Again, delayed SFHs give results very 
similar to those with declining SFHs. 

The predicted median SFR of all galaxies are compared in 
Fig- El As expected due to lower dust attenuation and older ages, 
the models with constant SFR show generally the lowest (cur- 
rent) star formation rates (cf. dBSS 12). As before, models A and 
E are very similar. Finally, assuming rising star formation histo- 
ries generally implies a higher SFR, due to, first, the lower Luv 
output per unit SFR (cf. Fig. 13, and, second, to a higher aver- 
age dust attenuation, needed to redden the UV slope, which is 
always dominated by the youngest stars. 



6 



D. Schaerer, S. de Barros, P. Sklias: Testing star formation histories and the SFR-mass relation of Lyman break galaxies 



-drops 



/ 

A 

, a* 



It**"* ** 1 
a* * " ■ 



WEI'.* 

pjjj^w^-* A A *^ " 4 



2 4 
XI - decl 



Fig. 8. Comparison of the median value of the reduced xl f° r 
each galaxy of the B-drop sample (z ~ 4) computed for different 
star formation histories and plotted as a function of xl obtained 
from model A (declining star formation histories). Same sym- 
bols and colors as in previous figures and listed in TableQ] 



3.2. Comparison of the fit quality 

Figure[8]shows the fit quality of the models expressed by the re- 
duced Xy value for the z ~ 4 sample. Overall, declining star 
formation histories (model A) give the best fit to most galaxies, 
and delayed SFHs very similar fits. Exponentially rising histo- 
ries provide better fits for some objects, but a higher^ 2 on aver- 
age. Models B and C with constant SFR and rapidly rising star 
formation yield fits of lower quality for most galaxies. 

We have also examined the influence of another reddenin g 
law (cf. Sect. El). App lying the SMC law of lPrevot et~ al. ( 1984); 
Bouc het et al.l ( 1 1985b we find that the majority of galaxies fit less 
well than with the Calzetti law, both the declining and rising star 
formation histories (model A, C, and D), although the differ- 
ences are generally relatively small (cf. Sect.|7]l. 



3.3. Implications for the SFR-mass relation and influence of 
the minimum age 

Figure|9]shows the SFR-mass relation for the 705 B-drop galax- 
ies as derived from the SED fits assuming the five different star 
formation histories. As already discussed in dBSS 12, the models 
show overall a large scatter and a m edian specific SFR ( SFR/M*) 
above the mean relation derived bv lDaddi et al.l (120071) at z ~ 2, 
also shown for comparison. In part the scatter in this diagram is 
determined by the possible range of SFR/M* predicted for stel- 
lar populations with different SFHs and different ages. Indeed, 
for a given L^i/Ly, probed observationally by Ljjy/Ly, the ratio 
SFR/M+ varies by ~ 2 orders of magnitude between young ages 
and t £ 1 Gyr for constant SFR models, and more for declining 



or delayed star formation histories (Figs. [2] and [3}. Of course, 
since reddening also modifies Lyv /Ly, this adds further scatter. 

As expected, the delayed SFH yields basically identical re- 
sults to the declining histories. Rising histories (model D) oc- 
cupy a similar domain as the other variable SFHs (models A, C, 
E), e xcept for the absen ce of galaxies below the mean relation 
from lDaddi et al.l (120071) . This is simply due to the fact that star 
formation cannot decrease in this model and for model B (con- 
stant SFR). Hence situations where the current SFR is low, but 
the UV flux is still present from somewhat older stars (e.g. B- 
type stars with lifetimes up to ~ 100 Myr) are not possible for 
these SFHs (cf. Fig. |2j. Objects well below the SFR-mass re- 
lation require obviously decreasing star formation histories, or 
some shut-off of star formation. 

The upper boundary of the SFR-M* plot, parallel to the 
indicated z ~ 2 relation, is determined by the minimum age, 
and by dust attenuation. Indeed, young ages imply lower masses 
and a higher SFR from simple consid erations of SFR/UV and 
mass/light ratios (cf. Figs. [2] and [3] also McL ure et al.ll201 ll) . 

Since age is kept free in all our model^l except for model B 
where we assume a minimum age of 50 Myr, and since some the 
median ages obtained are relatively young (e.g. compared to typ- 
ical dynamical timescales of disks; cf. dBSS 12) it is interesting 
to explore how the main physical parameters and the SFR-M* 
relation change when introducing a prior on the minimum age 
(or equivalently on the formation redshift) for all models. 

The effect of fixing the minimum age to — 50 Myr in the 
SED fits with exponentially declining and rising star formation 
histories (models A and D) on the SFR-mass relation is shown in 
Fig. [10] In short, for the declining SFHs the upper boundary be- 
comes now identical to the one obtained with SFR=const (model 
B), as expected, since declining SFHs include long timescales 
t for which the star formation history becomes undistinguish- 
able from constant, and since f m ; n = 50 Myr sets the maximum 
SFR/UV output. In contrast, for rising SFHs the upper boundary 
is found to be higher by a factor ~ 3 approximately. Again, this 
simply due to the less efficient UV output per unit SFR (or higher 
SFR/UV ratio) of the youngest stars, which always dominate the 
SFR for such star formation histories (Fig. [2j. The lower bound- 
ary of the SFR-mass relation remains of course unchanged by 
changes of the minimum age. In conclusion, introducing a mini- 
mum age affects in different ways the SFR-M* relation, depend- 
ing on the SFH. 

On average, i.e. for the entire sample of 705 B-drop galaxies, 
the median age and mass increase by a factor ~ 2, whereas SFR 
and Ay decrease by a factor ~ 2, when the prior f m j n = 50 Myr 
is introduced. For the rising histories (model D) the median age 
increases by a factor ~ 5, SFR decreases by 1.5, stellar mass 
increases by a factor 2. On the other hand the dust attenuation, 
both of individual galaxies and on average, remain very similar, 
with interesting implications on the predicted IR luminosity (see 
Sect.l4~4l 

To conclude this section we note that SED fits using models 
with variable star formation histories imply a fairly large scatter 
in the SFR-mass relation for LBGs at z > 3, as shown here and 
in dBSS 12. However, it must be pointed out that the amount of 
scatter cannot directly be translated into observable quantities 
such as UV and/or IR luminosity, which are commonly used as 
SFR indicators. Indeed, as we will show below (Sect. I4.31 >. a 
much smaller scatter is expected in terms of observed UV+IR 
luminosity. 



2 Here^J is defined as^ 2 /(A , Mt - 1), where 2Vgi t is the number of filters 
for which data is available for the galaxy, including non-detections. 



The only age limitation is given by the age of the Universe. 
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Fig. 9. Current SFR as a function of stellar mass of the z ~ 4 
LBGs obtained for different SF histories. Same symbols as Fig. 
|4] The dotted and dash-dotted lines delimit the area above which 
galaxies with Lir ^ 2. x 10 11 L Q (or equivalently a flux > 0.2 
mJy in ALMA band 7 ** which wavelength** for T c / = 35 K) or 
> 2. x 10 10 L Q , respectively, are found. The dashed line shows 
the SFR-mass relation at z ~ 2 for comparison. 



4. Predicted IR-mm emission 

4.1. IR luminosity and constraints on UV attenuation by dust 

Figure QT| shows the infrared luminosity predicted by the differ- 
ent models for the z ~ 4 sample, Fig.[T2]the corresponding his- 
togram of Lir/LuvU Since some galaxies have a median Ay = 
and hence no UV light is absorbed, their predicted median IR lu- 
minosity is Lir = 0. Such galaxies are not shown on these plots; 
their percentage ranges from 24% for model B (with constant 
SFR) to 10% for model C, which has the highest median Lir. 
Due to the adopted discretisation in Ay, the lowest, non-zero ex- 
tinction Ay = 0.05 corresponds to log(LiR/Luv) ~ -0.85, as 
seen e.g. in Fig. [26] 

As could be expected from the behavior of the UV attenu- 
ation, the models with constant SFR predict generally signifi- 
cantly lower Lir than models with different star formation his- 
tories. On_averag_e 1 _fhe highest Lir is predicted for steeply rising 
SF of iFinlator et alj d201 lh . followed by the exponentially ris- 
ing model D, and model A and E with very similar average IR 
luminosities. Clearly, if Lir can be determined with sufficient 
accuracy from observations, it should be possible to distinguish 
different models for individual galaxies, as well as statistically 
for large enough samples. 

Most importantly the IR/UV ratio allo ws one to mea- 
sure the UV attenuation a s well known (cf. iBuat et alJ [2005: 
Hi esias-Paramo et al.1l2007l and references therein), and this in- 
dependently of the assumed SF history and of the extinction law. 
To demonstrate this we compute the UV attenuation factor /uv 
defined by the ratio of the intrinsic (emitted) over the observed 



Fig. 10. Effect of adopting a minimum age prior on the SFR- 
mass relation. Green (red) show the models with declining (ris- 
ing) SFHs, as in Fig. [9] Including the prior t > 50 Myr yields the 
magenta (dark blue) points for the declining (rising) SFHs. Note 
the decrease of the upper boundary of the SFR-M* relation, but 
less so for rising SFH than for declining ones. Black triangles 
stand for model B with SFR=const. The dashed line shows the 
SFR-mass relation at z ~ 2 for comparison. 



UV luminosity, i.e. 



/uv - Luy/Lijv. 



(1) 



This quantity is shown in Fig. Q~3] for all galaxies of the z ~ 4 
sample, for all star formation histories and for the Calzetti at- 
tenuation and SMC extinction curve. The figure clearly shows a 
tight correlation between /uv and Lir/Luv independently of the 
assumed star formation history and of the ext inction/attenuatio n 
law, a behavior well known f r om the wor k of Bu at et al.1 J2001 : 
ll^lesias-Paramo et aD d2007l) : iBoquien et alJ d2012l) and others. 
The UV attenuation factor (defined here at 1800 A) can simply 
be described by 



log(/ uv ) = 0.24 + 0.44x + 0.1 6x 2 



(2) 



4 Here L uv is essentially model independent (except for the photo- 
metric redshift), since determined by the observed UV magnitude. 



where x = log(LiR/Luv)- For a specific attenuation law this can 
of course be translated into quantities such as Ay. For exam- 
ple, for the Calzetti law one has Ay = 2.5(Rylk^) log/uv = 
1.08 log/uv- Such a simple relation is quite naturally obtained, 
since the bulk of the energy reemitted in the IR after absorption 
by dust is also produced by the same stars contributing the bulk 
of the UV emission. More complicated relations could be ex- 
pected when e.g. older stellar populations such as AGB or RGB 
stars contribute t o dust emission, a s appears to be the case in low 
redshift galaxies dBuat et al.ll2005h . 

In short, measurements of Lir/Luv of individual LBGs pro- 
vide a direct measure of their UV attenuation, which should 
allow us to distinguish certain star formation histories, and in 
particular to test the systematic shift in physical properties pre- 
dicted between models assuming long star formation timescales 
(constant SFR) and variable (rising or declining) star formation 



8 



D. Schaerer, S. de Barros, P. Sklias: Testing star formation histories and the SFR-mass relation of Lyman break galaxies 



histories. However, additional observational constraints, such as 
emission line strengths discussed later, may be necessary to dis- 
tinguish between among variable SFHs, since both rising and de- 
clining histories yield fairly similar predictions for many physi- 
cal parameters (cf. dBSS12 and above). 
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Fig. 11. Same as Fig.|4]fbr the predicted IR luminosity. 



l°g(L, R Auv) 

Fig. 13. UV attenuation factor ftjy as a function of the predicted 
ratio of the IR/UV luminosity. A tight relation, fitted by Eq. Q] 
and shown by the dashed line, is obtained independently of the 
assumed SF history and of the extinction l aw. F or comparison 
the relation determined by iBoquien et alJ d2012h is shown by 
the dotted line. The corresponding attenuation Ay assuming the 
Calzetti law is shown on the right scale. 



-drops 



60 



40 



20 - 



H-lT 
Vv 



'-a 



EEL 



log(L IR /L uv ) 

Fig. 12. Histogram of the predicted median IR/UV luminosity 
ratio for the z ~ 4 sample from models assuming different star 
formation histories (Model A: green dashed line; model B: solid 
black; model C: solid blue; model D: red dashed; model E: yel- 
low dashed). 



4.2. IR and UV luminosity as tracers of the SFR in light of 
different star formation histories 

Often the IR luminosity or the sum of IR+U V luminosity are 
used as a measure of the star formation rate (e.g. 



Kennicuttl l998; 



Burgarella et all 120051; iReddv et all [20101 12013) 7 The basic as- 
sumptions made for such conversions are that of constant SFR 
and an age of ^ 50-100 Myr, necessary to reach an equilibrium 
output for the UV and bolometric luminosities (see also Fig. [2j». 
The determination of the SFR from Lir also assumes that the 
total (stellar) bolometric luminosity is reemitted in the infrared 
domain, which requires significant dust attenuation. Obviously 
the above assumptions cannot be valid in general, and in par- 
ticular they are not compatible with the different star formation 
histories explored in this paper. 

In Figs. [14] and [15] we show the predicted IR and IR+UV 
luminosity respectively as a function of the (current) SFR deter- 
mined from fitting our models to the sample of z ~ 4 galaxies . 
The "standard" SFR(IR) calibration given by iKennicutll fl998) 
- based on the bolometric luminosity - is also indicated for 
comparison. Overall quite large deviations from the simple cal- 
ibration are obtained for Lir — SFR, and smaller deviations for 
(Lir + Luv)-SFR. As expected model fits assuming constant SF 
(model B) agree very well with the Kennicutt relation when the 
IR+UV luminosity, providing a very good measure of the bolo- 
metric luminosity, is used. However, for galaxies with SFR ^ 
20-30 M yr , the predicted IR luminosity is decreased with 
the respect to this relation. This is simply due to the fact that 
these galaxies are not found to be dusty enough to reprocess a 
significant fraction of their UV light into IR radiation. This is 
also clearly illustrated by the large uncertainties in Lir, which 
can reach values down to if the pdf of Ay does so. This ef- 
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Fig. 14. Median IR luminosity Lir predicted for the B-drop sam- 
ple as a function of the median SFR obtained from the SED fits 
for the models with different star formation histories. Results 
from exponentially declining models (A) are plotted with green 
squares. The remaining symbols and colors are the same as in 
previous figures. For the constant SFR model the results, shown 
as black triangles, we also plot the 68% confidence range; for 
clarity these error bars are shown only for 1/6 randomly chosen 
galaxies of t he sample. The dotted line indicates the SFR(IR) 
calibration of lKennicutj d!998l) . 

feet, shown here for SFR=const, is present for all star formation 
histories. 

More generally, for varying star formation histories, both 
the IR and the IR+UV luminosity predicted by model fits are 
lower than expected from the Kennicutt relation between SFR 
and luminosity (Figs. [T4l and [15]). Two effects explain these dif- 
ferences. First, dust extinction needs to be high enough for Lir 
to be a good tracer of the total SFR (cf. above). Second, the ages 
are younger tha n the time neede d to reach the bolometric out- 
put assumed by iKennicutj (fl9 98). Typically, this is the case for 
ages f/r$ 1 for declining SFHs, somewhat longer for delayed 
and rising star formation histories, and -100 Myr for constant 
SF, as illustrated in Fig. [2] The second effect explains the val- 
ues of Lir + Luv below the canonical Lir-SFR relation. Added 
to that, the first effect further diminishes Lir for galaxies with 
lower SFR, as already explained above. In some cases, although 
for a small fraction of galaxies, we find that the predicted IR 
or IR+UV lum inosity is higher than what is expected from the 
iKennicuttl d 19981) relation using the model SFR. This is found 
with declining or delayed star formation histories for galaxies 
where f/r » 1, i.e. where star formation starts to decrease but 
longer-lived stars still contribute to the bolometric luminosity. 
These objects are expected to show weak emission lines, i.e. they 
should include galaxies from the category of "quiescent" LBGs 
identified by dBSS 12, and which can also be recognized by their 
3.6 micron excess, if at z ~ 3.8-5. 

Of course it should be noted that the SFR shown in these 
plots is not a direct observable, since it corresponds to the cur- 
rent SFR determined from the broad-band SED fits, and since 



Fig. 15. Same as Fig.[14]for the IR+UV luminosity, which pro- 
vides a better approximation of the bolometric luminosity. 



the SFR derived in this manner depends itself on model assump- 
tions, such as the assumed star formation history (cf. Sect.[3]and 
Fig-EJ- Observationally the best measure of the current SFR ob- 
tained from the SED fits should be the one measured by hydro- 
gen recombination lines (e.g. SFR(Hff)), once corrected for dust 
attenuation, since these emission lines trace best short-lived mas- 
sive stars. From the Ha flux predicted by our SED models, the 
consistency with SFR(Ha') can indeed by verified. 

The main conclusions to be drawn from Figs. [T4l and [T51 is 
that IR or (IR+UV) SFR calibrations assuming constant star for- 
mation over long timescales should be used with caution for the 
interpretation of high-z star forming galaxies (e.g. LBGs), since 
they are susceptible to variations of the star formation history 
on short timescales, for which there are good indications (see 
dBSS 12 and this paper). In any case, the effects discussed here 
need to be taken into account to properly predict the IR emission 
of LBGs and for comparisons with observations. Possible prob- 
lems related to the use of inconsistent SFR indicators, and im- 
plications on the SFR-mass relation are discussed below (Sects. 
l4~3landlT2i>. 

4.3. Implications on the SFR-mass relation 

Fig. [16] shows the SFR-mass relation with SFR transformed to 
direct observables, namely the sum of the UV and IR luminos- 
ity. Compared to the SFR-M* relation plotted in Fig. [9j which 
shows a very large scatter (with an upper bound up to 2 dex 
above the mean relation at z ~ 2) for models allowing for young 
ages and for variable star formation histories, the scatter in this 
relation is now significantly reduced, approximately by one or- 
der of magnitude, as judged by the upper boundary. The reason 
is again due to the fact that the models with the highest current 
SFR/M* values correspond to relatively young ages, where the 
bolometric output per SFR is lower (Fig. [2j. Hence the smaller 
dynamic range in bolometric luminosity, which is well approxi- 
mated by Lir + Luv- In other words, although the intrinsic, cur- 
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Fig. 16. Median UV+IR luminosity predicted consistently for 
the B-drop sample as a function of the median stellar mass ob- 
tained from the SED fits for the models with different star for- 
mation histories. Note that the predicted scatter in luminosity is 
much smaller than the one in the true, underlying SFR shown in 
Fig. [9] since the corresponding IR luminosity is computed con- 
sistently with the assumed SF timescales and ages. The dashed 
line indicates the expected UV+IR lumi nosity expected as sum- 
ing the SFR-mass relation at z ~ 2 from | Daddi et al.l d2007t) and 
the SFR(IR) calibration of lKennicuttl(ll998h . 



rent SFR may show a wide dispersion, the SFR one would infer 
from the UV+IR luminosity assuming "standard" SFR conver- 
sion factors would show a much smaller dispersion! This demon- 
strates that the finding of a small dispersion in the bolometric 
luminosity-M* plane does not exclude a high dispersion be- 
tween current SFR and stellar mass, and objects with very high 
values of the specific star formation rate, SFR/M* . 

4.4. Influence of the minimum age on the dust attenuation 
and IR luminosity 

How does the choice of a minimum age affect the derived atten- 
uation and hence the predicted IR luminosity? Figure [T71 shows 
the result of this exercise for the exponentially declining and ris- 
ing SFHs (model A and D). We find that introducing a lower 
limit ttmn (age prior) leaves the predicted IR luminosity of LBGs 
unchanged on average for the rising star formation histories, 
whereas Lir diminishes for declining SFHs. The basic reason for 
this behavior is the following: since for rising SFHs the intrinsic 
UV emission does not decrease when the age is increased, the at- 
tenuation and hence their IR luminosity remains comparable for 
most galaxies. In contrast, for decreasing star formation histories 
an older age allows one to reproduce the SED with a lower at- 
tenuation (cf. the well known age-extinction degeneracy). Thus 
a lower Lir, comparable to the values predicted for the constant 
SF model with the same minimum age, is expected for this case. 
In other words, although the predicted Lir depends in general 
on assumptions on the minimum age of stellar populations, this 
is not the case for rising star formation histories, which consis- 



Fig. 17. Effect of adopting an age prior (minimum age t > 50 
Myr) on the predicted IR luminosity for declining (green, model 
A) and rising (red, model D) histories. We plot the difference in 
log(LiR) between the unconstrained model and the one with age 
prior as a function of the latter. On average declining SFH show 
lower Lir whereas rising histories predict basically the same, or 
even higher IR luminosities. Rising SFHs case can thus be dis- 
tinguished from constant SFR, even if age priors apply. 



tently predict a higher UV attenuation and hence IR luminosities 
than constant SF or declining histories. A similar distinction was 
already seen in the SFR-mass relation (Fig. IT0T>. Measuring Lir 
can help test/constrain these different scenarios. 



5. Predictions for ALMA 

5.1. Continuum flux predictions 

In Fig. [18] we show predicted fluxes for the B-drop samples in 
one of the ALMA bands (band 7 covering ~ 0.8-1. 1mm) avail- 
able during the early ALMA cycles. As explained in Sect.[2]these 
fluxes are derived from the predicted IR emission, Lir discussed 
above, with the assumption of modified black bodies with a typ- 
ical dust temperature of Tj = 35 K and (3 — 2. Changes of the 
dust temperature by ±10 K correspond to a change of the flux 
by a factor ~ 2.5-2.8 approximately. At this redshift (z ~ 4) the 
emission increases (decreases) with decreasing (increasing) dust 
temperature in this band, since the peak of the spectrum lies at 
shorter wavelengths than band 7 considered here. 

Fig. [18] clearly shows the significant flux differences pre- 
dicted by models assuming different star formation histories. If 
the dust properties (temperature and exponent ft) do not change, 
the IR-mm fluxes simply reflect the differences in IR luminosity 
predicted by the models, as already discussed in Sect. |4] While 
the majority of the z ~ 4 LBGs show fluxes below the sensitivity 
limit reached with short (60 s) exposures with the full ALMA ar- 
ray, approximately half of the sample should be detectable with 
30min exposures. In any case, already with the current, incom- 
plete ALMA array the brighter part of the distribution should be 
detectable. Such measurements will provide badly needed, di- 
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Fig. 18. Predicted histogram of median dust continuum fluxes in 
the ALMA band 7 (0.8-1.1 mm) for the B-drop sample from 
models assuming different star formation histories. The predic- 
tions are plotted for modified black bodies with Tj = 35 K and 
P — 2.. Same color coding as in Fig. |2T] The vertical lines in- 
dicate the continuum sensitivity for an integration time of 60 s 
with the full array (long dashed), and with 32 antennae available 
for the ALMA cycle 1 (short dashed). The upper axis indicates 
the corresponding IR luminosity. Changes of the dust tempera- 
ture by +10 K correspond to a change of the flux by a factor ~ 
2.5-2.8 approximately. 



rect measurements of the UV attenuation of z ~ 4 LBGs and 
very useful tests/constraints on the SED models and on the star 
formation histories of these distant galaxies. 

When selecting galaxies by their IR luminosity, which part of 
the SFR-mass diagram such as Fig.|9]do we probe? Naively one 
would expect this to translate to a simple limit in SFR. However, 
since SFR/Lir is a priori not constant, i.e. age dependent, and 
since Lir also depends on the amount of dust attenuation, the 
limit is more complicated. This domain is e.g. illustrated for 
L IR > 2 x 10 11 L by the dashed line in Fig. [9] Whereas this 
criterium selects all massive galaxies with SFRS 35 M Q yr _1 , 
such galaxies with M 5= 6. x 10 8 M have relatively young popu- 
lations and less dust, which implies a lower IR luminosity for the 
same SFR, below the IR luminosity limit. Although the detailed 
shape of curves of LiR=const in Fig.|9]depend somewhat on the 
value of Lir, their qualitative behaviour is similar. 




Fig. 19. Predicted IR luminosity as a function of the absolute UV 
magnitude for the B-drop sample (z ~ 4) assuming declining star 
formation histories (model A). Colors show the full 2D probabil- 
ity distribution function (pdf) of the sample. The solid (dashed) 
lines encircle the 60 (90) % confidence range. The median val- 
ues of each galaxy with a median A y > 0.05 are shown as white 
circles. The median of the pdf in bins of UV magnitude is shown 
as red stars. The dotted line shows, for comparison, the expected 
Lir wh en assuming the standard SFR(UV) and SFR(IR) rela- 
tions of lKennicutdfi998h . and SFR(UV)=SFR(IR). 
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Fig. 20. Same as Fig. [19] for the i-drop sample with 60 galaxies 
at higher redshift (z ~ 6). Note the significant decrease in the 
median Lir-Muv relation compared to z ~ 4, due to the decrease 
of the median dust attenuation found with increasing redshift (cf. 
dBSS12). 



5.2. Evolution of the IR luminosity with redshift 

The typical predicted IR luminosities at z ~ 4 are shown in Fig. 
[19] as a function of the absolute UV magnitude (here M1500 for 
comparison with dBSS 12) for model A, i.e. declining star forma- 
tion histories. The figure shows the 2D probability distribution 
function (pdf), the median values for each galaxy (white circles), 
and the median of the pdf in bins of UV magnitude (red stars). 
As mentioned earlier, the pdf also includes a tail to very low IR 
luminosities, set somewhat arbitrarily to Lir values outside the 



plot, since some of the SED fits also allow for zero extinction. 
This explains why the median IR luminosity (red stars) is shifted 
towards lower values than the peak of the pdf. Incidentally, at 
z ~ 4, the median Lir for model A is found to be quite close to 
the IR luminosity predicted from the UV m agnitude using th e 
standard SFR(UV) and SFR(IR) relations of iKennicuttl 0998), 
and assuming SFR(UV)=SFR(IR), as shown by the dotted line. 
On average we find that the UV brightest galaxies are also ex- 
pected to be brightest in the I R, in agreement with iReddv et al.l 
(l2012h : lGonzalez et alJd2012l) and other studies. Typical median 
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luminosities are of the order of Lir ~ 10 ll) 11 L Q for model A 
(cf.Fig.HS). 
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Fig. 21. Histogram of the predicted median (rest-frame) Ha 
equivalent widths for the z ~ 4 sample from models assum- 
ing different star formation histories (Model A: green solid line; 
model B: solid black; model C: solid blue; model D: red dashed; 
model E: yellow dashed). 

In Fig. [20] we show the IR luminosity for the i-drop sam- 
ple (z ~ 6) predicted from model A. It should be noted that a 
higher fraction of the 60 objects of this sample are found with a 
median attenuation of zero (below our minimum non-zero value 
of Ay = 0.05), implying thus a median Lir outside this plot. In 
other words the complete pdf has a secondary peak at low values 
of Lir, not shown here, which explains - as for Fig. [19]- the shift 
of the median of the sample (red stars) with respect to the appar- 
ent peak of the pdf. Clearly, the predicted IR luminosity of z ~ 6 
LBGs is lower by a factor ~ 3-A at most UV magnitudes, due 
to the lower dust attenuation at high redshift (see dBSS12 for 
more details). However, despite the lower reddening, the SED 
fits show that galaxies at z > 6 with dust exist quite likely, even 
at these high redshifts (cf. dBSS12, also Schaerer & de Barros 
2010). Indeed, these studies show a clear trend of increasing av- 
erage dust attenuation with galaxy mass, at all redshifts. 



6. Predicted (rest-frame) optical emission lines 

Since one of the main features of our SED analysis is the treat- 
ment of nebular lines, whose effects turns out to be not negli- 
gible for approximately two thirds of the LBGs (cf. dBSS12), 
it is obviously of interest to examine the predicted strength of 
these emission lines, and to test these predictions observation- 
ally. From the 63 lines from the (rest-frame) UV, optical and 
near-IR domain included in our models, we illustrate here Ha, 
among the strongest optical lines, and [O n] /13727, the bluest 
optical line. Their predicted equivalent width distribution for the 
B-drop sample and for the different star formation histories are 
shown in Figs. [2T| and [22] 



250 



200 



150 




100 



log(EW([0II]3727)) 

Fig. 22. Same as Fig.[2T]for the (rest-frame) [O n] A3121 equiv- 
alent width. 
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Fig. 23. Predicted (rest frame) Ha equivalent width of the B-drop 
galaxies as a function of the stellar mass. Different colors indi- 
cate the different models, as in previous plots. For all SFHs an 
anti-correlation of Wn a with galaxy mass is expected. However, 
the detailed values of the equivalent width of each galaxy, and 
the range of Wh u covered, are model dependent. 



The predicted Ha equivalent widths cover a wide range, with 
a maximum reaching up to ~ 3000 A for some rare cases of very 
young galaxies. Most SFHs, except for the model with constant 
SFR, predict a fairly wide distribution of equivalent widths. The 
rising star formation histories (models C and D) predict higher 
median equivalent widths, since the current star formation al- 
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Fig. 24. Predicted correlation between the (rest-frame) equiva- 
lent width of Ha (large symbols) and the specific star formation 
rate. The values shown here are derived for the sample of B-drop 
galaxies. However, all samples and models yield very similar re- 
lations. Small symbols show the same for EW([0 n] /I3727), 
which traces less well the sSFR. 

ways produces Lyman continuum photons, hence also emission 
lines. As expected, models with constant SFR predict a rela- 
tively narrow range of equivalent widths, limited by a minimum 
Whq.^ 100 A, which is reached close to the maximum age of $ 1 
Gyr. Only models allowing for some decline of star formation 
(e.g. model A, E), can predict objects with very weak or no emis- 
sion lines, as clearly shown in Fig.l2"Tlandl2"2"l From these figures 
is evident that measurements of the equivalent width distribution 
in LBG samples can, in principle, distinguish different star for- 
mation histories. In practice, however, such line measurements 
are currently difficult or not yet possible, due to limited sensi- 
tivity/access at the observed near-IR wavelengths. For this task, 
recombination lines such as Ha should ideally be used, as they 
trace directly current star formation, whereas metal lines depend 
more strongly on nebular conditions (metallicity etc.). 

The Ha fluxes predicted forz ~ 4 sample are between 10" 18 
and 10~ 16 erg s _1 cm 2 , with the lowest median predicted for 
the SFR=const model. The predicted [O n] A3721 fluxes are in a 
similar range, but less differences are expected between different 
SFHs, as already mentioned above. 

Not surprisingly, we find a relatively strong (anti-)correlation 
of the Ha equivalent width with galaxy mass, shown in Fig. [23] 
Obviously different models occupy somewhat different regions 
of this plot, but in all cases do we find such a behavior. Indeed, 
this trend is expected, since the average age decreases and the 
specific star formation rate (sSFR=SFR/M*) increases with de- 
creasing galaxy mass for all SFHs (dBSS12). Wn a decreases 
also with galaxy age, but the spread in age can be very large, 
given various star formation timescales. The best correlation is 
indeed found between the Ha equivalent width and the specific 
SFR, shown in Fig. [24] Again, this is a natural trend, since the 
Ha equivalent width measures the ratio between the Ha flux — 
tracking the current SFR — and the underlying continuum flux, 



which depends on the accumulated stellar mass. We note that all 
models predict a very similar relation, which also holds for all 
samples from z ~ 3 to 6, although they occupy different regions 
of the relation. Measurements of Who (from spectroscopy or es- 
timates from broad-band photometry) may therefore be used to 
determine the sSFR of individual high-z galaxies. As Fig.l24lalso 
shows, metal lines such as [O n] /13727, which are not recombi- 
nation lines, trace the sSFR less well. 

Our predicted equivalent widths are comparable to those 
found recently in emission line galaxies at z ~ 0.6 - 2.4 for 
[O ill Xilll, fa ml /U4959,5007, and Ha dStraughn et al.ll2009l: 
lAtek etalJbOl U Ivan der Wei et al.ll201 ll) An anti-correlation of 
Wu_ a with galaxy mass is found at all redshifts, where data is cur- 
rently available, from z — to 2-2.6, as recently discussed by 
Fumagalli et al] d2012l) . Furthermore, the predicted Ha equiva- 
lent widths of the B-drops are i n good agr e emen t with the obser- 
vation s of z ~ 2-2.6 LBGs bv lErb etal] d2006l) : iMancini et all 
d201lh in the same mass range (typically at M* 3; 10 10 M Q ). Our 
models thus predict that this trend of increasing Ha equivalent 
width continues down to LBGs of lower mass. 

The 47 strong Ha emitters identified by I Shim et aT 1 (120111)1) 
at z ~ 4 are included in our LBG sample, and amply discussed in 
dBSS12. Their selection, based on an excess in the 3.6 /vm filter 
with respect to 4.5 yum, is sensitive to EW(Ha)^ 200 A. Clearly, 
our method is able to detect/predict weaker lines, although of 
course with increasing uncertainty (not shown here). T he high- 
est inferred equivalent widths from Sh im et alj d201 lbl) . well in 
excess of 1000 A, are comparable to our values, as are the pre- 
dicted Ha fluxes. Although very rare at low redshift, galaxies 
with EW(Ha) > 500 A also exist in nearby galaxies, as pointed 
out by[Sh lm &Charvld2012h . They are all found at masses well 
below log Mi, S= 9.5 M and at high sSFR, quite compatible with 
our model predictions. 

We conclude that the emission line strengths predicted by 
our models from fits to the observed broad-band photometry and 
the correlations found are compatible with the currently existing 
data and trends. Future observations of emission lines in z ~ 3-6 
LBGs with existing and new ground-based facilities, and with 
the James Webb Space Telescope will hopefully provide further 
tests on the importance of emission lines and constraints on the 
star formation histories and specific star formation rate of these 
galaxies. 

7. Discussion 

7. 1 . Influence of different extinction laws 

All our SED fits have been carried out using the Calzetti attenua- 
tion law, commonly thought to be appropriate for starburst galax- 
ies. To examine the effect of other laws we have rerun all our fits 
for the B-drop sample with the ste e per attenuation / extinc tion law 
for the SMC iPrevot et all dl984l) : iBouchet et all dl985l) . Other 
laws, intermediate between the relatively grey Calzetti law and 
the steep SMC curve, should yield intermediate results. 

Overall, we find a lower fit quality (higher reduced xl) f° r 
the vast majority of objects using the SMC law. However, for 
most of them the difference remains relatively small (A%1 ^ 1). 
We verified that no systematic effect is found as a function of 
galaxy age, as could e.g. be expected if young galaxies (e.g. t < 
100 Myr, as judged by the age derived assuming the Calzetti 
law) would be better fit w ith the SMC law, as suggested e.g. by 
iReddv et al.l(l2010il2012l) . 

Since the SMC law is steeper than that of Calzetti, SED 
fits for LBGs are generally yield a lower UV attenuation /uv 
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Fig. 25. Results/predictions from SED fits of the B-drop sample 
for model A (declining SFHs) assuming the Calzetti attenuation 
law (green squares) or the SMC extinction law (magenta). Top: 
IR+UV luminosity as a function of the stellar mass. Bottom: 
Current star formation rate versus stellar mass. 



(and even more so a lower Ay), older ag es, and hence on av- 
erage somewhat larger stellar masses ("e.g.lPapovich et"ai1 l200 1 ; 
IVerma et al J 120071: lYabe et al.ll2009l iReddv et al.ll2010l) . As ex- 
pected, our fits including nebular lines also yield the same av- 
erage trends. The young ages favored with the Calzetti law for 
some objects by the strong Ha excess at 3.6 pm are in part com- 
pensated by increasing the star formation timescale t for declin- 
ing SFHs. The resulting Ha equivalent width distribution (cf. 
Fig. [2Tb remains quite similar, with the largest differences due to 
the assumed star formation histories. 

Important differences of the SED fits with the SMC law are 
in particular a lower SFR and less IR emission. The former is 
due both to the older age (hence lower SFR/UV, cf. Fig. |2]l and 
to the lower UV attenuation. The lower attenuation implies im- 
mediately a lower IR luminosity. For model A (declining SFHs) 
the median SFR of our B-drop sample is lowered by a factor 
~ 2, from 11.6 to 5.4 M yr _1 , the median IR luminosity from 
4. x 10 10 L by a factor ~ 5, compared to SED fits with the 
Calzetti law. The effect of adopting the SMC law on the SFR- 
mass on the (LiR+Luv)-mass relations for model A are shown in 
Fig. [25] The top panel clearly shows the reduced IR luminosity, 
implying a lack of B-drop galaxies with Lir ^ 2. x 10 11 L Q (in 
this sample) if the SMC law was applicable. A large scatter is 
still found in the SFR-mass relation with some objects reaching 
large specific SFRs, as shown in the bottom panel. However, at 
masses ^ 10 9 M the SED fits with the SMC law predict ba- 
sically no LBGs with high sp ecific SFR, i.e. above the z ~ 2 
relation of lDaddi etaD (120071) . This shift to lower SFR values 
is again mostly due to age shifts, leading to ages ^100 Myr 
for all galaxies with M 10 9 M . It must, however, be remem- 
bered that the model SFR shown here represents the current SFR, 
which shows a wider dynamic range (spread) than the observable 
UV and IR luminosity (cf. Sect.l431l. 



Fig. 26. Predicted IR/UV luminosity ratio as a function of the ob- 
served UV slope determined from SED fits of the B-drop sample 
for model A (declining SFHs) assuming the Calzetti attenuation 
law (green squares) or the SMC extinction law (magenta). Black 
triangles show the predictions for model B (constant SFR) with 
the Calzetti law for comparison. The solid (dashed) lines show 
the expected relation for galaxies with an intrinsic UV slope of 
B k -2.2 for the two attenuation laws. 



At low redshift, and ocasionally at z ~ 2-3, observations of 
the IR/UV ratio and measurements of the UV slope B have been 
used to distinguish the attenuation law (e .g. iBaker et afl feOOl: 
ISiana et alJl2008ll2009tlReddv et alj|2010h . The so-called IRX- 
B plot, based on the predicted IR luminosity plus the observed 
Luv and B, for our z ~ 4 LBG sample and for the two atten- 
uation laws is shown in Fig. [26] As already mentioned above, 
a lower IR/UV ratio is found from SED fits using the SMC law. 
Besides an accurate measurement of the IR/UV ratio (soon feasi- 
ble with ALMA), this method requires also a good knowledge of 
the intrinsic UV slope, and accurate measurements of B to con- 
strain the attenuation law. Indeed, as clear from this plot, UV 
slope measurements of individual galaxies show a large scat- 
ter and their errors are relatively large since the B slope de- 
pends on one UV col or, here B oc 5.3 x (i - z) color, following 
Bou wens et al. (2009). More accurate measurements are feasi- 
ble, e.g. using more photometric bands or ba sed on SED fits 
dFinkelstein et : al.ll201 it ICastellano et al.ll2012l) . In any case, the 
use of the IRX-B plot to distinguish attenuation laws also re- 
quires the knowledge of the intrinsic UV slope, which may not 
be unique, since a priori dependent on the star formation history, 
age, and metallicity of the stellar population. 

Distinguishing different attenuation laws may thus not be 
straightforward. However, some combinations of star formation 
history and attenuation law are clearly distinguishable once IR 
luminosities are measured with ALMA. As already mentioned, 
our SED fits to the present data yield in general better fits with 
the Calzetti law for the vast majority of objects, although the 
difference in^ is not very large. 
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7.2. Comparison with earlier studies 

We shall now briefly discuss other papers, which have presented 
tests for or constraints on different star formation histories for 
Lyman break galaxies using also the IR as a constraint. 

To the best of our knowledge, our study is the first discussing 
SED fits with various star formation histories including both de- 
clining and rising parametrisations for a large number o f z ~ 
3-6 LBGs. At somewhat lower redshift Reddy et al. (2012) have 
recently examined a sample of ~ 300 LBGs with MIPS 24 fim 
detections and spectroscopic redshifts 1.5 < z < 2.6. The 24 fim 
flux is used to infer the IR luminosity, which in turn provides 
an independent constraint on the SFR and dust content (as dis- 
cussed abo ve and in ear l ier pa pers, cf. lReddv et ail 2008). From 
their work iReddv et al.l (1201 2l) conclude that rising star forma- 
tion histories provide a better agreement between different SFR 
indicators, and hence that declining star formation histories may 
not be accurate for typical galaxies at z ^ 2. 

Obviously, this conclusion is in contrast to ours, since we 
find that declining SFHs provide the best fits for the majority 
of z ~ 3-6 LBGs, and that delayed and rising star formation 
histories fare comparably. How can we understand these differ- 
ences? In their comparison of star formation rates Redd v~et~aT] 
(120 12h find two different behaviors. First, for galaxies classi- 
fied as "young" (f < 100 Myr) according to their SED fits with 
the Calzetti law, they find that the SFR inferred from the fits, 
SFR(SED), is systematically larger than the SFR(IR+UV) de- 
rived fro m the UV a nd IR luminosity using the standard calibra- 
tions of iKennicutd (Il998l) . From this they conclude that there 
is an apparent conflict with young ages, which they propose 
to resolve by invoking a different, steeper attenuation law for 
young objects. We question the solidity of such a conflict, and 
argue that their resolution is by no means unique. Indeed, al- 
though they note that a young age is not compatible with as- 
sumptions made for the standard UV-SFR conversion, they ne- 
glect the fact that both IR and UV conversions depend on the 
SFH and age in a similar way (see their Fig. 25 and Fig. [2]), 
modifying the (UV+IR) — SFR conversion upward by typically 
twice the amount they state. The effect they find is clearly shown 
to predicted for our z ~ 4 sample, as illustrated in Fig.[l5j where 
we note that for the majority of objects fitted with unconstrained 
ages and found to be younger than 100 Myr from our models, the 
SFR(SED) is larger than the SFR(IR+UV) obtained from an in- 
consistent application of the standard calibrations. A consistent 
analysis of the time and star formation history dependent UV 
and IR emission may th erefore also s i mply explain the apparent 
problem pointed out by Red dv etaf] ([2012) for young galaxies 
without resort to an age dependent attenuation law. 

The second different behavior noted bv IReddv et all (1201 2h 
is for galaxies classified by their model as being older than 100 
Myr. For these they find that exponentially declining SFHs yield 
SFR(SED) systematically lower than SFR(IR+UV), from which 
they conclude that such histories are disfavored. Indeed, their 
models apparently yield best-fit solutions with f/r » 1 for these 
galaxies, which nat urally explains the discrepancy between these 
SFR indicators (cf. lReddv etaIll2012L and Fig.©. However, UV 
continuum observations, showing clear signatures of young mas- 
sive stars, and Ha m easurements for som e of the LBGs indicate 
f/r < 1 according to lReddv et al.l d2012l) . From this, the authors 
conclude that exponentially declining SFHs are statistically in- 
co nsistent with the observations. We cannot verify the results 
of IReddv et al.l (120121) as their observational data is not pub- 
lished. However, we note that age is the most model dependent 
quantity and that systematically younger ages are obtained when 



nebular emission is taken into account (cf. above and dBSS12). 
Furthermore in our analysis of z ~ 3-6 LBGs we find values 
of f/r < 1 for the majority of object s, in agreement wit h the 
observational arguments mentioned bv lReddv et al.l (120121) . We 
may thus suspect that a reanalysis of z ~ 2 LBGs with models 
including nebular emission may yield different results. 

In short, we find that the conclusion of IReddv et all (I20I2I) 
favoring clearly rising star formation histories over declining 
ones for z ~ 2 LBGs is not robust. In any case, stellar popu- 
lations/SFHs of LBGs at higher reds hift may be different, and 
from the result of IReddv et al.l (1201 2l) at z ~ 1. 6-2 A one cannot 
conclude that exponentially declining SF histories are not appro- 
priate for LBGs at z > 3. It is therefore important to obtain direct 
measurements of the IR luminosity at higher redshifts. 

Other comparisons with Re ddv et alj d2012l) regarding con- 
stant and rising star formation histories have already been dis- 
cussed in dBSS12, and shall not be repeated here. 

7.3. Possible caveats 

Before closing this section, we wish to remind the reader that 
our SED fitting tool suffers from the same possible limitations 
as other models, including e.g. the unknown star formation his- 
tories of galaxies (modeled here with simple parametrisation), 
the question of multiple populations and the "outshining" prob- 
lem due to young populations, uncertainties related to the ex- 
tinction (which law? a unique attenuation for all stars? etc.), un- 
certainties in the stellar initial mass function, and of course also 
uncertainties in stellar evolution and at mosphere model s . Som e 
of these issues ar e e.g. addressed by iMaraston et al.l d2010l); 
IWuvts et alj d2009l): ICharlot & Longhettil (|2001|) ; IWalcher et all 
(1201 lh : lLevesque et al.l (120121) : ISchaererl d2012l) . 

In this respect our results should of course be considered 
as "differential", i.e. with respect to other models making the 
same/similar assumptions, but neglecting the effects of nebular 
emission and of variable star formation histories, which are ex- 
plored here systematically for the first time for a large sample of 
Lyman break galaxies. Concerning nebular emission our models 
make some simple assumptions, which may also influence the 
results. Our main assumptions are no loss of ionizing photons 
(i.e. no Lyman continuum escape, no dust inside the Hn regions 
etc.), and the same attenuation for the stellar continuum and neb- 
ular lines. Both assumptions maximize to some extent the possi- 
ble contribution of nebular emission, and could therefore lead to 
an overestimate of the effects of nebular emission on the derived 
physical parameters. The results would then be intermediate be- 
tween the present case with nebular emission and the models 
without, discussed in depth in dBSS12. Exploring variations of 
the above assumptions may be interesting and worth pursuing in 
future studies. In any case, it is clear that for a majority (approx- 
imately two thirds) of LBGs the inclusion of nebular effects im- 
proves the quality of the SED fits, and that models without neb- 
ular lines cannot reproduce the observed excess in certain photo- 
metric bands (most prominently the 3.6 yum excess for galaxies 
with redshift between 3.8 and 5). For the bulk of the LBGs our 
models should therefore provide a better description than previ- 
ous models neglecting nebular emission, despite the simplifying 
assumptions made. To progress further on these questions addi- 
tional observational constraints, such as direct measurements of 
IR luminosities, medium-or narrow-band photometry, or mea- 
surements of emission lines should be of considerable help. The 
predictions presented here should also provide a useful base for 
future comparisons and tests of the SED models. 
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8. Summary and conclusions 

Following up on our earlier detailed study (Ide Barros et al .2012) 
of a large sample of LBGs from redshift z ~ 3 to 6 located in 
the GOODS-South field, using for the first time an SED fitting 
tool including the effects of nebular emission on the synthetic 
photometry, we have examined the impact of different star for- 
mation histories (SFHs) on the derived physical parameters of 
these galaxies, on the SFR-mass relation, on different SFR in- 
dicators (UV and IR), on the expected dust extinction and the 
corresponding IR luminosity, and on the expected strengths of 
emission lines such as Ha and [O n] A3727. 

To do so, we have carried out SED fits for five different SFHs 
including exponentially rising and so-called delayed SFHs, plus 
the three histories already considered in dBSS12 (see Table Q] 
and Fig.Q}. Metallicity is also treated as a free parameter, and we 
have examined the effect of two different extinction/attenuation 
laws (Calzetti and SMC). The usual physical parameters, stel- 
lar mass, SFR, age, and attenuation are derived from the SED 
fits to the broad band photometry reaching from the U band 
to 8 fan, using Monte Carlo simulations to derive their median 
values (and the detailed probability distribution function, gener- 
ally not discussed here). We have also computed consistently 
the predicted IR luminosities, Lir, for all galaxies, assuming 
energy-conservation, i.e. that all the radiation absorbed by dust 
is reemitted in the IR. Finally, the predicted IR luminosities have 
been translated to flux predictions in various IR bands, assum- 
ing modified black body spectra. The Lir predictions allow us in 
particular to examine in a consistent way the effects of variable 
SFHs and ages on this observable quantity, showing thus signif- 
icant departures from results assuming inconsistent SFR(IR) or 
SFR(UV) calibrations. 

Our main results concerning the impact of star formation his- 
tories on the physical parameters of LBGs, exemplified to a sam- 
ple of 705 LBGs at z ~ 4 (B-drop galaxies), can be summarized 
as follows (see Sect. [3}: 

- Compared to commonly adopted SED fits assuming con- 
stant SFR, no nebular emission and an age prior of t > 50 
Myr, models with exponentially declining SFHs, nebular 
lines and no age constraint yield younger ages, lower stellar 
masses, higher current SFR, higher specific star formation 
rates (sSFR=SF R/M»), and higher dus t extinction (Ay), as 
already shown in de Barr os etaD d2012l) . 

- Exponentially declining SFHs yield overall the best fits for 
the majority of LBGs. 

- Assuming delayed star formation histories one obtains basi- 
cally identical physical parameters (and fit qualities) as for 
exponentially declining SFHs. 

- Rising star formation histories with variable timescales im- 
ply generally a similar stellar masses, and comparable or 
somewhat higher dust extinction than models assuming de- 
clining SFHs. The latter leads to the highest star formation 
rates and to similar or higher IR luminosities as for declining 
histories. 

- Overall "standard" models assuming constant and neglecting 
lines predict systematically higher stellar masses, lower ex- 
tinction, lower SFR, lower IR luminosities, and more narrow 
range of equivalent widths for optical emission lines than all 
the other star formation histories considered here. 

Combining these different physical parameters we obtain the 
following (Sect. |U: 

- We find significant deviations between the derived SFR 
and IR luminosity from the commonly used SFR(IR) or 



SFR(IR+UV) calibration of iKennicuttl (Il998l) . Such differ- 
ences naturally arise, due to differences in the derived ages 
and in the adopted star formation histories. In most cases (i.e. 
for most galaxies and SFHs) we find that the Kennicutt re- 
lation will underestimate the true, current SFR derived from 
the SED fits (Sect. 14.21 1. Consistent SED studies including 
also the IR are therefore necessary, if the SFHs and ages may 
differ from those assumed in standard SFR calibrations. 

- A large scatter is found in the SFR-mass relation for mod- 
els with declining an d delayed SFH and no age prior (cf. 
Ide Barros et alj|2012l) . The same also hold for models with 
rising star formation histories. The scatter is reduced when a 
minimum age (e.g. t > 50 Myr) is adopted. Even in this case, 
models with rising SFHs can show a large scatter, since high 
SFRs are found to the high(er) extinction. 

- A large scatter in the SFR-mass relation does not necessarily 
imply the same scatter in the Lir (or (Lir + Luv)-mass rela- 
tion and vice versa, when the IR luminosity is computed con- 
sistently from the chosen SED model (i.e. accounting for age 
and SFH effects). The same also applies to SFR(UV)-mass 
diagrams where the UV luminosity and the corresponding 
standard SFR conversion is used. We suggest that the true 
scatter in the SFR-mass relation obtained in this way may 
indeed be underestimated, if the true star formation histo- 
ries are variable on relatively short timescales. Indeed such 
SFH variations can reproduce more successfully features re- 
lated to emission lines, such as the observed 3.6 /mi ex- 
cess in z ~ 4-5 LBGs. They may also be more relevant to 
at higher redshift, where the dynamical timescales decrease 
with(l +z)- 3/2 . 

Our consistent predictions of IR luminosities (and fluxes) 
show that different SFHs lead to significantly different amounts 
of reddening and hence to different IR/UV luminosity ratios. 
Measurements of IR luminosities of individual LBGs or statisti- 
cal samples of such galaxies can be used to distinguish different 
SFHs, and hence also different specific SFRs predicted by such 
models. ALMA observations will thus be able to provide inde- 
pendent constraints on behavior of the sSFR at high redshift, and 
on the scatter in the SFR-mass relation. 

We show predictions for the IR luminosities of B-drop and i- 
drop galaxies for different star formation histories and as a func- 
tion of UV magnitude. The typical/median Lir is predicted to 
be Lir ~ 10 10 " 11 L for LBGs with absolute UV magnitudes of 
M uv ~ -22 to -19. 

Finally we also show the predicted strengths (equivalent 
widths) of the Ha and [O n] /13727 emission lines. Again, dif- 
ferent star formation histories naturally lead to different EW dis- 
tributions, which can in principle be used to constrain the SFHs. 
Our models predict on average higher equivalent widths in low 
mass galaxies, in agreement with currently available observa- 
tions at z < 3, and a clear anti-correlation of Who with the spe- 
cific SFR. 

Our predictions should in particular provide new tests us- 
ing IR observations with ALMA and/or measurements of (rest- 
frame) optical emission lines to obtain a better insight on the star 
formation histories of high redshift LBGs, on the behaviour of 
the SFR-mass relations and on the evolution of the specific SFR 
with redshift. 
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